Nonlinear Spectral Singularities and Lasing Threshold Condition 
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A spectral singularity is a mathematical notion with an intriguing physical realization in terms 
of certain zero-width resonances. In optics it manifests as lasing at the threshold gain. We explore 
the application of their recently-developed nonlinear generalization in the study of the effect of 
a nonlinearity on the lasing threshold condition for an infinite planar slab of gain medium. In 
particular, for a Kerr nonlinearity, we derive an explicit expression for the intensity of the emitted 
waves from the slab and discuss the implications of our results for the time-reversed system that 
acts as a coherent perfect absorber. 
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I. INTRODUCTION 

One of the most interesting properties of complex scat- 
tering potentials [l| is that unlike their real analogs they 
can admit scattering states satisfying Siegert outgoing 
boundary conditions [|[ . Because these states have a real 
and positive energy, they behave exactly like a zero- width 
resonance. This observation is originally made in [H in 
an attempt to understand the physical meaning of the 
mathematical notion of a spectral singularity [|[ . It has 
motivated the study of physical aspects and applications 
of spectral singularities particularly in optics |3|, [5j-tl0| . 

A remarkable outcome of this study is that for an op- 
tical system consisting of an infinite planar slab of ho- 
mogeneous gain medium the well-known lasing threshold 
condition follows from the equation ensuring the gener- 
ation of a spectral singularity @. Furthermore, if we 
consider the time-reversed setup, where the slab consists 
of a lossy medium with the loss coefficient having the 
same value as the gain coefficient required for a spectral 
singularity, it behaves as a coherent perfect absorber Q. 

A more recent development is a nonlinear general- 
ization of the notion of a spectral singularity [11] . A 
characteristic feature of nonlinear spectral singularities 
is that they are intensity-dependent. This means that 
they correspond to the emission of waves with a particu- 
lar wavelength-amplitude profile. Ref. [ll| offers details 
of this phenomenon for a delta-function potential with a 
complex coupling constant 

In order to improve our understanding of the properties 
and possible applications of nonlinear spectral singular- 
ities, in the present article, we study their behavior for 
an infinite planar slab of gain material subject to a weak 
nonlinearity. In particular, we give a detailed analytic 
description of the effects of a Kerr nonlinearity. 

Consider a linearly polarized electromagnetic wave 
that is aligned along the x-axis in some cartesian co- 
ordinate system and has an electric field of the form 
E = e~ tuJt £(z)e x , where u> is the angular frequency of 
the wave, $ is a complex-valued function, and e x is the 
unit vector along the positive x-axis. Suppose that this 
wave interacts with an infinite planar slab of gain mate- 



rial of thickness a that is aligned in the x-y plane and 
located between the planes z = and z = a. Then it is 
well-known that the Helmholtz equation for this system 
reduces to the Schorodingcr equation, 

-S"'(z) +V{z)£{z) = k 2 S{z), (1) 

for the energy-dependent complex barrier potential [§| 

V(z):=k 2 (l-n 2 ) X a(z), (2) 

where k := lo/c is the wavemmiber, n is the complex 
refractive index of the slab, and 



Xa(z) := 



1 for < z < a, 
otherwise. 



(3) 



When a nonlinearity is present, n 2 is replaced by n 2 + 
gf(\S'(z)\) and ([1} takes the form 

-S"( z ) + V{z)S(z)-gk 2 Xa{z)f(\S{z)\)S{z) = k 2 S{z), 

(4) 

where g is a real constant and / is a real- valued function 
1131. In terms of the scaled variables, 



z/a, 



ip(x) := $ (ax), M.:=ak, 

j2/i „2\ 



(5) 



x 

7 :=-il 2 g, 3:=£ 2 (l-n 2 ), »(x) := 3Xi(x) 
we can write (QJ as 

-V"(x)+t;(x)V(x)+7Xi(x)/(|V(x)|)V(x) = & 2 yj(x). (6) 
For a Kerr nonlinearity, where 

/(H) = h/f, (7) 

([6]) is known as the Gross-Pitacvskii equation and has 
applications in the description of the Bose-Einstein con- 
densates [3 • For real values of 3 it admits exact analytic 
solutions in terms of the Jacobi elliptic functions [15| ■ 

II. NONLINEAR SPECTRAL SINGULARITIES 

In Ref. [ri| we introduce nonlinear spectral singular- 
ities associated with In short, they correspond to 



2 



the real and positive values of .ft 2 for which ((6j) admits a 
nonzero solution fulfilling the outgoing boundary condi- 
tions: 

V>'(0) + i£V(0) =0, ^'(l)-i^ fc (l) = 0. (8) 

In order to see the physical implications of these con- 
ditions, we use the fact that for x ^ [0, 1] the general 
solution of © is a linear combination of the plane waves 
e ±lfoc . This allows us to envisage scattering solutions of 
© corresponding to incident waves from the left (x < 0) 
and the right (x > 1). The former has the form [ill ] : 

N_ e -i$* + N_ e ix* f or x < o, 
ip l (x) = I C(x) for < x < 1, (9) 

N+e iMx for x > 1, 

where £ is a solution of © that satisfies 

C(l) = N + e m , C'(l) = iW + e m , (10) 

N + is a complex constant that gives the amplitude of the 
transmitted wave, and 



iGJjZ) 
2ft 



iV_ 



iG+{&) 
2ft 1 



G±(&) :=C'(0)±^C(0). 



(11) 
(12) 



We can use © and (fTTj) to derive the following expres- 
sions for the reflection and transmission coefficients from 
the left, respectively fill ]. 



R l = - 



G-(S) 



2ikN+ 



(13) 



In view of (fT0|) the solution (j9)) satisfies the second equa- 
tion in ([8]). The first of these equations, which is equiv- 
alent to G + (K) = 0, implies that both the reflection and 
transmission coefficients diverge. This is a characteristic 
property of a resonance state. Notice however that the 
energy ft 2 does not have an imaginary part. Therefore 
© is a scattering solution of (|6]) which behaves like a 
resonance with a zero width Q. 

Now, suppose that we tune the parameters of the 
system so that it supports a spectral singularity with 
wavenumber k and amplitude parameter N + . Then any 
left-incident plane wave with wavenumber and ampli- 
tude, respectively, in a close vicinity of k and AL. is am- 
plified and emitted from both sides of the slab. In the ab- 
sence of the nonlincarity, this phenomenon is amplitude- 
independent and the system amplifies the background 
noise and serves as a laser that functions at the very 
threshold gain Q . The presence of the nonlinearity mod- 
ifies the lasing threshold condition by introducing in it a 
particular amplitude-dependence. Our main objective is 
to derive an explicit expression for this modified lasing 
threshold condition. 



III. LASING THRESHOLD CONDITION 

Determination of nonlinear spectral singularities re- 
quires the solution of the boundary-value problem de- 
fined by ([5]) and (JSJ). This is equivalent to finding the 
function £, which satisfies © and ([TU]) . and demanding 
that 



G+(.ft) =0. 



(14) 



For x G [0,1], we can express ^ in the form of the 
integral equation 

V>(x) = Vo(x) + 7 f G ( x - y)/(l^(y)l)^(y)rfy, (is) 



where t/'o is the general solution of the linear equation 



ip" + n 2 & 2 i> = 0, 



(16) 



in [0,1], x £ [0,1] is arbitrary, and G is the Green's 
function for this equation, i.e., 



G{u) := 



sin(n^w) 



(17) 



Repeated use of (fl"5j) in its right-hand side yields a per- 
turbative series expansion for ip with 7 playing the role 
of the perturbation parameter. In particular, we find the 
following first-order perturbativc expression for 



C(x)«Co(x)+ 7 Ci(x), 



(18) 



where sa means that we neglect quadratic and higher or- 
der terms in powers of 7, £0 is the solution of (|16p satis- 
fying (TIUj) . and 

Ci(x):= f G(x- y )/(|Co(y)|)Co(y)dy. (19) 



It is easy to show that 



Co(x) 



N+e iM 



2n 



(20) 
(21) 



(n + l) e ^(*-i) + ( n _ !) e -*nJi(x-l) 

Next, we insert (fT5|) in (TT2"j) to obtain 

G±(J?)«Gi 0) (n,il)+7Gi 1) (n,^), 
where 

G ( £(n,R) := Cj-(0) ± i^0(0) for j = 0, 1. (22) 
The computation of G± ^ (ft) is straightforward. It gives 



G (0) M = Ar +e '-*ft(n 2 -l)sin(nft) ^ 



(23) 



Gf(n,ft) 



2n 



L(n,&), (24) 
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where 



where 



L(n,R) 



-2inM 



(25) 



According to (f2~4")l and (|25l) , in the absence of nonlinearity, 
7 = 0, and we find a spectral singularity provided that 
the right-hand side of (|25|) vanishes. This is in complete 
agreement with the results of 

Let rto and be such that L(no,^o) = 0, i.e., they 
correspond to a linear spectral singularity, and 770 and 
kq be respectively the real and imaginary parts of no, so 
that no := ?/o + * K o- Then in view of (|25p and the fact 
that for typical gain media, 



r) - 1 > -/co > 0, 



we have 



-ino-So 



1p - 1 
n + 1 ' 



(26) 



(27) 



Recalling that the gain coefficient go necessary for the 
emergence of the spectral singularity is related to kq ac- 
cording to go = — 2^oKo/a, and taking the modulus- 
square of both sides of (|2"T)) . we obtain 



,9o 



2a ln |ft| 2 ' 



(28) 



where Tl := (n - l)/(n + 1). In view of ([H, K ~ 
(i]o ~ l)/( 7 ?o + 1)- Therefore, it gives the reflexivity of the 
slab, and (f2"5|) coincides with the standard expression for 
the lasing threshold condition Q, 



go « - In - 

a \t]o - 1 



(29) 



In order to obtain nonlinear corrections to ([2U)) , we 
solve ([14]) perturbatively. Let ni and £1 be such that, up 
to linear terms in 7, Eq. (fl4| holds for 



n +7"i> 



fL= &o +7^1. 



(30) 



Inserting Q30p in (|21[) , expanding the resulting expression 
for G + in powers of 7, and equating the coefficient of 7 
to zero give 

9 no Gf (n ,^o)ni + ^ Gf (no.^o^i+G^Cno,^)) = 0. 

(31) 

Eq. ([27| simplifies the calculation of d no G® (no, So) and 
9i? G^(no,^o) enormously. The result is 



n 



(33) 



The calculation of (no, So) is more involved. How- 



ever, it turns out that we can use (H 
write it in the form: 



, and (J27J) to 



G^(n ,Mo) = -N + e^cl [ f (\N + c no h(x)\)h{xfdx, 

Jo 

(34) 



n + l 
2n 



h(x) 



-irio.Hoa: 



For example, consider a Kerr nonlinearity ([TJ. Then we 
can easily perform the integral in ([34|) and use (|27|) to 



express (tio,So) as 



G?'(noA) = 



8i\N + \ 2 N+e liio {Anl 



3) 



^(9n4 + nS 4 -10|n | 4 ) 



(35) 



We can determine the nonlinear corrections to the lo- 
cation of the spectral singularities by substituting (|32[) . 
(f3"3"|). and ([55)1 in (|3"Tj) and solving the resulting equation 
for any two of Re(ni), Im(ni), and Si in terms of the 
third. 

In order to gain a better understanding of the conse- 
quences of ([3"Tjl. we expand the right-hand sides of (|32j) . 
(|3"3"|) , and (f55|) in powers of «o and keep the leading order 
term. This gives 



9 no Gf (n ,^o) 

0*, Gj (no, So) 
G^K,^) 



A^°Sq [(7?g - 1)S) - 2i] 

m 

N+e m °MVo - 1), 
3\N+\ 2 N + e iRo (ril - 1) 
4776 s -So K o 



, (36) 
(37) 
(38) 



The gain coefficient g that is required to support a 
nonlinear spectral singularity is given by 



2 &K 



9o 



1 + 7 h -- 

K -KQ 



(39) 



where n :— Im(n) and ki :— Im(ni). Suppose that we 
are interested in a nonlinear spectral singularity corre- 
sponding to a Kerr nonlinearity that has the same wave- 
length as its linear analog. Then &i = and we can 
use (|3"Tj) and (|3"6"|) - (|3T))) to determine the modified lasing 
threshold condition: 



Ho 



9o 



1 - 



6(772 - 1) 7 |JV 4 



6g|7V+| 2 



4] ln 2 ( 



2( r)o + l 
)7o-l 



77 2 (77o 2 -l)ln 2 (^±i) 



(40) 



Here in the last relation we have switched to the unsealed 
nonlinearity parameter g and use the fact that for typical 
optical setups (770 — l) 2 /So ~ 0. 

Eq. (|4"0)) is our main result. It implies that for a posi- 
tive Kerr nonlinearity parameter g, which is usually the 
case [IH, the presence of the nonlinearity increases the 
necessary gain for producing a spectral singularity. More 
importantly it shows that the slab acts as an amplifier for 
a left-incident plane waves only if g exceeds go- In this 
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case it emits an amplified transmitted plane wave whose 
intensity is given by 



9o 



12g5o 



VoiVo - l)ln 



( m + 


t)] 


\Vo - 





(41) 



Because this relation follows from our first-order pertur- 
bative analysis, it applies for situations where g|A + | 2 <C 
1. Typically g < 1CT 13 cm 2 /W which for a value of | 1 2 
as large as 1 GW/cm 2 gives g|iV + | 2 < 1CP 4 . This shows 
that (j41j) is quite reliable. 

We can view (|4T|) mathematical justification for 
the well-known fact that lasers operate when the gain 
coefficient is larger than go- It is important to notice 
that in the linear case the spectral singularity disappears 
as the value of the gain coefficient exceeds go and the slab 
stops functioning as a laser. This is not the case when we 
take into account the effect of a Kerr nonlincarity with 
g>0. 

We have performed a first-order perturbative calcula- 
tion of G-(&) and used the result together with (fTTj) 
to compute the amplitude parameter 7V_ for the emitted 
wave from the left-hand side of the slab. This calculation 
shows that whenever we arrange the parameters of the 
system to generate a nonlinear spectral singularity, the 
contribution of the terms of the order 7 to N_ cancel one 
another and we find A_ ps N + e lix . Therefore similarly to 
the linear case, the wave emitted from the left-hand side 
of the slab has the same amplitude and phase as the one 
emitted from its right-hand side. For the time-reversed 
system, this implies that a lossy slab serves as a coher- 
ent perfect absorber provided that the loss factor a has 
a numerical value that is larger than go and the intensity 



1 7V_)_ 1 2 / 2 of the incident coherent waves be given by (f4Tj) 
with g replaced with a. 



IV. CONCLUDING REMARKS 

An interesting outcome of the search for physical impli- 
cations of spectral singularities is a mathematical deriva- 
tion of the lasing threshold condition. This is quite differ- 
ent from the standard derivation of this condition which 
is essentially based on the physical principle of conserva- 
tion of energy (l6j . In the present article we have used 
the newly developed concept of a nonlinear spectral sin- 
gularity to account for the contribution of nonlinearities 
to the lasing threshold condition. 

Our method is quite general and can be applied 
to almost all nonlinearities of physical interest. For 
a Kerr nonlincarity it provides a simple explanation 
for the well-known fact that our slab model docs 
not function as a laser at the threshold gain go- It 
starts emitting radiation only when the gain coeffi- 
cient g exceeds go- In this case, the intensity of the 
emitted radiation turns out to be proportional to g — go ■ 
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